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Abstract 

We prove in this paper that the electro-weak SU{4:) l ® U{l)x gauge models 
with spontaneous symmetry breaking can offer a natural framework for generating 
neutrino masses by simply exploiting the tree level realization of dimension-five 
effective operators. The novelty of our approach resides in the fact that the scalar 
sector needs not to be enlarged, since these operators are constructed as direct 
products among scalar multiplets already existing in the model. There is a unique 
generic matrix for Youkawa couplings in the neutrino sector. The charged leptons 
are already in their diagonal basis. This framework can lead to a suitable fit of the 
established phenomenology for the left-handed neutrinos, while the right-handed 
neutrino masses come out in the sub-fceV region, independently of the cut-off A. 
The latter introduces in the theory an intermediate scale (however, more close to 
GUT than to SM) at about lO^^GeV which is a crucial ingredient for the left- 
handed neutrino phenomenology. 

PACS numbers: 14.60.St; 14.60.Pq; 12.60.Fr; 12.60.Cn. 

Key words: dimension-five effective operators, see-saw mechanism, 3-4-1 gauge 
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1 Introduction 

One of the main challenges in particle physics 1 1 1 today is the neutrino mass issue. Both 
its origin and order of magnitude are still awaiting a compelling theoretical explanation. 
Observational collaborations |2| such as SuperKamiokande, K2K, SNO, KamLAND, 
LSND and others have stated the phenomenon of neutrino oscillation as an indisputable 
evidence. This state of affairs claims for tiny but non-zero masses for neutrinos, re- 
gardless they will prove themselves as Dirac or Majorana particles. Consequently, the 
theory is called to supply a convenient framework and suitable mechanisms for gener- 
ating these tiny masses, assuming that the Standard Model (SM) does not necessarily 
include right-handed neutrinos (otherwise unavoidable ingredients in accomplishing a 
non-zero mass term). The experimental side of the neutrino mass issue enforces certain 
restrictions, namely the observed mass splitting ratio ta ~ Ar7iQ/ATO^(„j ~ 0.033 
and particular patterns for the mixing angles Oq ~ 34° and 6q ~ 45°, along with a 
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likely 613 ~ 0. The absolute mass hierarchy remains still undetermined on theoretical 
grounds. What we only know at present is that it lies in the sub-ey region| 1 1. 

We mention the see-saw mechanism|3 | (with its variants) and higher-dimension 
operators in effective theories [4] among the most appealing theoretical devices de- 
signed to accommodate neutrino phenomenology and predict viable consequences of 
it at low energies. These approaches generally require a larger framework than the one 
offered by the SM. For instance, the canonical seesaw mechanism essentially relies on 
a new higher scale (not subject in the gauge symmetry of the SM, but violating B — L 
symmetry) to generate Majorana right-handed masses. In a self-explanatory notation, 
the seesaw 6x6 matrix can be put as: 

/ m]^ 

MM+d ^ 

\ TTlD Mr 

By diagonalizing it one can get both the left-handed and right-handed neutrino 

masses as 

A'l (lyL) = -mlM^^mo , M (lyn) ^ Mr. (2) 

Therefore it seems worthwhile to investigate some extensions of the SM that in- 
clude in a natural way right-handed neutrinos. During the last two decades, gauge 
models such as SU{3)c «) SU{3) l ® U{l)x f^"^ ) and SU{i)c ® SU{^)l® U{l)x 
ITJ - L9J have been intensely investigated. They offer promising results when address- 
ing the neutrino issue. In most of the 3-3-1 models Q a scalar sextet must be added to 
the Higgs sector in order to give rise to some seesaw mechanisms. In 3-4-1 models the 
same strategy was considered in Refs. |8| where a new Higgs decaplet is introduced 
in the scalar sector of the model in order to set up proper seesaw terms. Yet, this can 
affect the boson masses previously calculated, while some ad hoc hypothesis on the 
breaking scales of the decaplet must be speculated. 

Our approach here avoids such new ingredients and proves itself able to set up the 
canonical seesaw mechanism via dimension-five effective operators constructed out of 
the existing ingredients in the model. It deals with the class of 3-4-1 models without 
exotic electric charges, so that even the exotic new quarks exhibit the ±| and ±i 
electric charges and mix with the traditional quarks of the SM. It also does not involve 
any enlargement of the scalar sector, but just makes use of the cut-off A up to which the 
model works as a safe renormalizable effective theory. The latter not only determines 
the highest bound for the validity of the theory but also play a crucial role in predictions 
regarding the order of magnitude of the active neutrinos' masses. 

Our paper is conceived as follows: after a brief overview of the gauge model and 
its main phenomenological features in Sec. 2, we focus on its lepton content in Sec. 
3 presenting the mass generating procedure based on dimension-5 effective operators 
in the neutrino sector. Sec. 4 is reserved for some numerical estimates and predic- 
tions in a particular scenarios taken into consideration, while in Sec. 5 we sketch our 
conclusions. 
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2 Electro-weak SU{i) l ® U{l)x gauge model 



We start here by briefly presenting the particular 3-4-1 gauge model under considera- 
tion here. However, the reader can find its phenomenological details treated in extenso 
in Refs. |8|. Evidently, our focus will go to its lepton sector, as we intend to exploit 
the realization of dimension-five effective operators responsible for giving rise (at the 
tree level) to the well known seesaw terms. 

In the electro-weak gauge group (5[/(3)c ® SU{4)l ® U {l)x the electric charge 
operator is a linear combination of diagonal Hermitian generators from Cartan sub- 
algebra. It is realized in the manner: Q = T^l + -^Tf^L + -^^isl + \XI where 

Ta = \\a normalized as Tr{TaTb) = ^Sat- 
The lepton representations in this model are: 
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For the sake of completeness, we present, in addition, the quark representations 
of the 3-4-1 model of interest here. In order to cancel all the chiral anomalies, two 
left-handed quark families must transform as Qa — { D'^ )^ ^ 

(3,4, —1/6), differently from a third generation which does it in the manner Q3 = 

(da U3 U U' )^ ~ (3, 4*, 5/6). Their right-handed partners are singlets with 
respect to the electro-weak gauge group, namely dsij, daR, D^r, D'^ri.^ (3, 1, —2/3), 
U3R, UaR, Ur, U'jj^ (3, 1, 4/3), with a = 1, 2. Capital letters denote exotic quarks, 
yet their electric charges exhibit (in the particular class of models under consideration 
here) the same pattern as ordinary quarks. Note that other classes of 3-4-1 models Q 
allow for new quarks with some exotic electric charges such as ±4/3 or ±5/3. How- 
ever, the exact colored symmetry SU{3)c of the QCD remains to describe the strong 
interaction as a vector-Uke theory, and its predictions are not affected at low energies 
by this enlargement of content. 

The gauge bosons of the electro-weak sector occur in connection with the standard 
generators TaL of the su(4) algebra. In this basis, the gauge fields are of U{l)x 
and S su(4), that is 
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(5) 



withDi^ = Al/V2 + Al/V6 + AiyVT2, D^^ = -AI/V2 + AI/V6 + A^^Vu, 
= -2AI/V6 + Ajj^/Vu, Di^, = -iAflVvl as diagonal Hermitian bosons. 



3 



By inspecting Eq.©, one notes that - apart from the charged Weinberg bosons (VT^) 
- there are several heavy degrees of freedom, namely two new charged bosons 
and along with and (and their complex conjugated). The diagonal 

entries provide us with the neutral physical bosons: the massless photon A™' for the 
electromagnetic interaction and massive Z^, along with two heavy neutral bosons Z'^ 
and Z'' involved in the neutral currents of the model. 

The scalar quadruplets in order to break the symmetry of the model stand in the 
following representations: 



Xi_ 

V ) 

( \ 

V j 



(1,4,-3/2) 



\ a ) 



(1,4,1/2) 



(6) 



(7) 



The superscripts denote their electric charge in units of e and i = 2,3,4. 

The electro-weak sector of the model provides us with two distinct couplings: g 
for SU{4:)l and gx for J7(l)xrespectively. Hence, the covariant derivatives read: 
Df_, = d^- i{gA^ + gx§A°^). Evidently, g is the SM coupling of the SU{2) l and g' 
of the C/(l)y, since SU{2)l ® U{1)y must be a subgroup of 5C/ (4) l ® C/(l)x- With 



P 
P\ 



respect to this subgroup, a Higgs doublet occurs from (f>^ ' - namely p - 
(1,2,1/2) - and two Higgs singlets p^, (1,1, — 1/2)- Consequently, there is 
also ^ ~ ( - 1 ^ (1; 2, 1/2) and corresponding singlets X2 ' X3 ^ (1, 1, ^V^)- 
When the symmetry is spontaneously broken up to the SM electro-weak group, these 



two scalar doublets can be seen as the traditional 



P\ 



(1,2,1) of the 



SM - with the completely decoupled p^, p^^ (1, 2, 0) at this level - and (p = ia 
(1, 2, —1), due to the equivalence 2 ~ 2* specific to 5C/(2) only. 



X 
Xl 

Now, the symmetry breaking pattern becomes quite obvious. The four scalar mul- 
tiplets in Eqs. (I6]l - dTji break the symmetry of the model in three steps to the residual 
one, namely to the electromagnetic U{l)em- 



SU{3)c(E> SU{4)l(E>U{1) 



X 



SU{3)c <E) SU{3) L <E) U{1)^' 



V 



V + v'^ 



SU{'i)c® SU{2)l®U{1)y 
5C/(3)c ;7(l)em 



by developing the vacuum expectation values (vev): < 0^^^ >= ( u' )^ 
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, < 0(2) >= ( V )'^, < 0(3) >^ ( V < 0W >^ 

(O y)^.A reasonable alignment is assumed here t; = w' ^ y = y 
in order to get rapid and simpler estimations. Evidently, v, v' — 174GeV is the 
SM electro-weak breaking scale, while the new scales V, V are specific to this 3-4-1 
model. 

In the symmetry breaking limit one can easily obtain the couplings match, namely: 

11 1 

This relation leads straightforwardly to the 

g% ^ sin^ 9w{mx) 

When applying the renormalization group procedure to Eq. dSOl l in the limit sin^ 6]^ 
sin^ OwiiTiz) it results that the unification scale mx is not sensitive to ax{mz) which 
remains ax{mz) ^ 10^^ even for m^up to lO^^GeV (while a{mz) ^ 1/128) . 

Of special interest is the boson mass spectrum. It can be computed as: 



m 



^Ki^lg'iV^ + v'^) , fn%o^KO-.)^l9Hv' + v^) , (10) 
mjc± ^ ^g^r' + v'^) , mj,o^xo*)^lg^V'' + v^) , (11) 



cos-^ Ow 2 



M'z'Z" - TTffM • (12) 



In our approximation V = V , one can easily diagonalize the above matrix and 
get: m\, = \g^V^ and m|„ = g^V^. Z^'' couples only with exotic fermions having 
no interaction with ordinary particles and thus totally decoupling from the low energy 
phenomenology, while Z' mixes with the SM Z and exhibits a mass (17] [D) lower 
bounded by 2TeV (or greater), in close dependence on which one of the three gener- 
ations of quarks transforms differently from the other two. The restrictions on the Z' 
neutral boson come from plugging into this model the data supplied by atomic parity 
violation experiments and certain processes in the meson systems 13° ~ D^, — B^, 
B^ — B^, — (with their corresponding mixings in the FCNC). However, these 
theoretical bounds Q [8] are consistent with recent experimental observations |[T] that 
suggest a lower bound around ITeV for Z'. 

In order to ensure the consistency with the SM phenomenology, the heavy particles 
of the 3-4-1 symmetry must not compromise the precision tests of the SM, particu- 
larly the oblique corrections must remain unaffected. That is, the new heavy quarks. 
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new heavy bosons and the extra scalar fields must give negligible contributions to the 
oblique parameters U, S, T, namely to the vacuum polarization amplitudes of the W 
and Z bosons. Being completely decoupled from low energy physics, Z'^ does not 
contribute to these oblique corrections. Also, yO(yO*) does not contribute as long it 
is a singlet under SU{2). The SU{2) doublets {K+,K°)and iX+,X°) (their com- 
plex conjugates) does not alter the one-loop calculations due to their degenerate masses 
(see Eqs. ( fTOt - (fTTTi). There remains to be investigated only the contribution of some 

scalar doublets such as I I I I whose couplings are subject to a proper 

tuning. However, these detailed calculations exceeds the aim of this paper 

All the fermion masses are generated through Higgs mechanism by scalar particles 
interacting with fermion fields. In this connection, the large split of breaking scales 
works efficiently in preserving the SM phenomenology. Ordinary quarks and leptons 
acquire their masses at the SM scale v while exotic quarks at the new V, V scales. 

A vast amount of theoretical research has been accomplished in this field due to 
several striking features that make such SM-extensions quite appealing and much valu- 
able despite the fact that they claim for a plethora of new particles. We count some of 
their assets, (i) First of all, the generation number in the fermion sector seems to get its 
explanation (absent in the SM, where simply an ad hoc triplication of the first genera- 
tion is performed). In order to cancel the chiral anomalies - this time by an interplay 
among families - the number of generations must be divisible by the number of colors 
Nc = 3. This leads to exactly 3 generations if one assumes the asymptotic freedom 
condition from QCD that limits them to no more than 5. (ii) Contrary to SM, these 
models supply a natural framework for charge quantisation (Doff and Pisano, |7 1). (iii) 
The strong CP problem can be elegantly solved due to a natural existence of Peccei- 
Quinn symmetry (Pal, Montero, Sanchez-Vega |5|). (iv) It offers a suitable framework 
for implementing the little Higgs mechanism (Kong [7] ). (v) All SM neutral currents 
and masses are identically recovered, (vi) All new particles acquire their masses from 
the high scales V, V so they do not interfere with the SM phenomenology at low en- 
ergies supplied by present facilities, (vii) If the third generation of quarks is the one 
transforming differently, that accounts naturally for the unbalancing heaviness of the 
top quark . 

We mention that our work focuses on a particular model from the 3-4-1 class of 
models, namely the one corresponding in the systematic classification accomplished 
by Ponce and Sanchez |i9J to 6 = c = 1, Model A. Its rich phenomenology of such 
models can be found and compared in Refs. QE). 



3 Lepton masses 

The Yukawa sector of any gauge model is set up to supply fermion masses (conse- 
quently the SSB). There are introduced certain Yukawa coefficients (matrix h) that 
couple left-handed and right-handed fermion fields. We write down the most general 
combinations allowed by the gauge symmetry. 
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Z:^^* = h\J^L<P^^^ + kR + (h^J S+f^L + hp+f^^ + hp'+f<r^) + h.C. (13) 

where = ei,/ii,TL. S* matrices are defined as follows Sr — {4>^'^^ ® tp'^^^ + 
<^(4) ® ^ (1, 10, 1), Sd = (<^(2) ® (/'(^^ + 0(3) - (1, 10, 1), S'o = 

(0(2) ^ ^(4) _|_ 0(4) ^ 0(2)>) ^ ^j^^ j^Q^ 2.) with h^j as thegeneric Yukawa matrices 
for Dirac and Majorana terms. Evidently, this is the basis where the charged leptons are 
already diagonal, so h\j — 0. Up to this point we proved that the electro-weak gauge 
symmetry SU{4)l <S) U{i)x allows for a natural implementation of neutrino masses 
via dimension-five effective operators plus the Yukawa couplings, usually subject to 
certain extra assumptions to overcome their arbitrariness. 

In order to restrict ourselves to a simpler version the entries in the generic neutrino 
matrix are the same regardless their nature, namely hfj = hf^ Therefore, the Yukawa 
Lagrangian becomes in the case at hand here: 



h.C. (14) 



In concrete expressions below the Yukawa coefficients will be denoted in order 

-A = hee, B = h^^, C = hrr, D — hep,, D' = h^e E = her, E' = hre, F = h^r, 
F' = hrp. 

It is natural to consider that the positions 3 and 4 in each lepton quadruplet are 
precisely N = vr and N' = i/^. As long as they are sterile with respect to Z and 
exhibit indistinguishable couplings to the new Z' and Z" bosons (see for instance Ref. 
ISJ) This assumption leads straightforwardly to the following identification: 



— ^Y 

By inspecting Eq.(fT3Tl one can easily identify 



rD' 
L,Y . 



(15) 



m(e) = /igU, to(/z) 



TO(r) = /i^w . 



(16) 



since we work in a basis where the charged lepton sector is diagonal. 

Taking into consideration the field theory mass formulas for respectively Dirac 
(£y — —■m-o'^^i\) + i/.c.) and Majorana {Ly = —\niM^'^'i> + i?.c.) terms, along 
with the above vev alignment and the normal seesaw (|TJ, one gets the mass matrices: 



Mu = h—, 



(17) 



T/2 



Mr= (h + h^)—. 



(18) 



From Eq. (|2]l, in the flavor basis the Majorana terms for left-handed and right- 
handed neutrinos can be read 



M{vl)^ h^{h + h^) h 



V 



(19) 
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(20) 



By diagonalizing Eq. ( ITsT i with a proper Ur one obtains those matrices in the mass 
basis. First step is: 



{h + hT) = Diag (ri, r2, r3) ^Ul{h + h^) Ur 
Now, inserting this resuh in Eq. (fT9] l one computes 



h^UR{h + K^) ^Ul^h 



V 

T 



which can be diagonalized as Ml — U'^ M{vl)U , namely: 



Ml = 



U^h^UR{h + hT) ^UlhU 



V 

T 



(21) 



(22) 



(23) 



The physical neutrino masses can be computed via Eq.(|23]l if we consider their mix- 
ing (for details, see the reviews in Ref. ifTOl ). The unitary mixing matrix U (U^U = 1) 
links the gauge-flavor basis to the physical basis of massive neutrinos: 



(24) 



i=i 



where a = e, (corresponding to neutrino gauge eigenstates), and i = 1,2,3 (corre- 
sponding to massive physical neutrinos with masses nii). The mixing matrix U pmns 
(Pontecorvo-Maki-Nakagawa-Sakata) has in the standard parametrization the form: 



(C12C13 S12C13 sise" 

-S12C23 - Ci2Sl3S23e** C12C23 - Sl2Sl3S23e*'' C13S23 
S12S23 - Ci2Si3C23e*'' -C12S23 - Sl2Sl3C23e** C13C23 

(25) 

where the notations sin(?23 = S23, sin0i2 = S12, sin0i3 = S13, cos023 — C23, 
cos 012 — C12, COS013 — Ci3 Stand for the mixing angles and S is the Dirac CP phase 
(with phenomenological meaning). Usually, to this one a diagonal Majorana phase 
matrix P = Diag {l,e^°',e^^) is sticked, though it can be absorbed by redefining 
fields. The standard identification leads to solar angle - 612, atmospheric angle - ^23, 
reactor angle - 613. Several different patterns for matrix|25]has been considered in the 
literature. The most appealing approach stemming from the seminal work of Harrison 
Perkins, Scott ifTTl . the so called "tri-bi-maximal" ansatz is largely invoked when the 
PMNS matrix is analysed and its phenomenological consequences are worked out in 
different models. A possible alternative to it followed the "bi-maximal" line and was 
developed in Refs. |[T2| . These particular textures - in good agreement with data - are 
often unfolded by enforcing certain discrete flavor symmetries on M. 

The global data 1 1 1 regarding neutrino oscillations impose certain restrictions, namely 
sin^ 9i2 ~ 0.3, sin^ 6*23 — 0.5, and a small reactor angle 6'i3(probably near zero), along 
with the mass sphttings Atoi2 — 7.6 x lO^'^eV"^ and Aml^ ~ 2.4 x 10~^eV'^. 



8 



4 Plausible scenario and numerical estimates 



In this section we prove that our construction is not a mere theoretical device. On 
contrary it can work very well when it comes to confronting the experimental data. We 
don't claim to make a general analysis here and get general predictions, but just apply 
the above procedure to a very particular scenario to get a quick fit with the data. The 
results make it obvious that our approach can be further developed and employed in 
investigating the neutrino sector. 



4.1 Left-handed physical neutrinos 

Since Ur is not restricted on observational ground, one can assume for the sake of 
simplicity a suitable scenario in which hji = —hij - namely D' = —D, E' = —E, 
F' = —F - that leads to Ur = I, and hence: 



{h + h^)=Diag{2A,2B,2C), {h + h^) ' = Diag{^, ^) (26) 



and consequently: 
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Bearing in mind that Trace is independent of the basis we work in, so that TrM{vi,) = 
m^L, one obtains: 

If all the Yulawa couplings are in the same range with the coupling of the charged 
rlepton (at most comparable but no greater than it) then one obtains an upper bound 
for the sum of the individual neutrino masses. This is: 

Y,miL<\m{T)^ (29) 

i 

As long as the experimental evidence imposes an upper bound on the range of left- 
handed neutrino masses (a few eV), one can estimate the cit-off energies up to which 
this model is valid. lfv= 174Gey and m(T) = 1777Mey, then A > lAxW^GeV 
which evidently is an intermediate level between SM and GUT energies. 
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4.2 Mass hierarchy and mixing angles 

In order to get some rough estimates of the mass spectrum and fit properly the mixings 
between the three species of left-handed neutrinos, let's consider a simple setting: A ~ 
-0.371 X 10"^ B ~ C ~ -0.5383 x 10~^ D E ^ 0.9034i x lO'^, F ~ 
— 0.17266i X 10~^, where i"^ ^ —\ the complex unity. 

Under these circumstances, the left-handed neutrino mass matrix (IZTT i becomes: 

/ 0.150 0.145 -0.145 \ 2 
M{vl) ^ 0.145 1.350 1.100 x 10"^^ (30) 
\ -0.145 1.100 1.350 / ^ 

so that it can be roughly diagonaUzed by the unitary matrix 

/ 0.831 0.556 \ 
U= -0.393 0.587 0.707 (31) 
\ 0.393 -0.587 0.707 / 

corresponding to the experimentally observed values sin^ 612 ~ 0.31, sin^ 6*23 — 0.5, 
andsin^ 0i3 ~ 0. 

The mass spectrum comes out in this particular case in a normal hierarchy, namely 

|toi| ~ 0.00216ey, ~ 0.08866eF, ~ 0.52983eF. (32) 

The mass splitting ratio yields ta — 0.03 in good agreement with experimental 
data and so the mass squared splittings are. 

Of course, a further work could consider a non-zero CP phase violation and a small 
but non-zero reactor angle so that from this stage one can perform a more accurate 
calculus and take into consideration a plethora of scenarios once the solar (^12) and 
atmospheric angles (6*23) are firmly established. All these scenarios, of course, cax 
lead to many different M{vl) matrices and, finally, even to a different mass hierarchy. 
Our results here are nothing but the proof that our method can work and is not claimed 
to be a general analytical analysis (which will be performed in a future work). 

4.3 Right-handed sterile neutrinos 

Now, a few words about the right-handed neutrinos in our approach framework. They 
acquire (Eq.(l26ll) the following masses in the scenario at hand: 

91/2 

= — (A + B + C) (33) 

■i 

According to the assumption in Sec. 4.1, this can be approximated as: 

> ^ ™« - 6m(T)— (34) 

i 

or equivalently: 
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^ (7) (35) 

which, as expected, is not affected by the cut-off A. If V is not very high - say around 
1 ~\QTeV - so that its new physics is testable at LHC, the so called sterile neutrinos de- 
velop masses in the suh-keV range. The stability of these three species of right-handed 
neutrinos could recommend they as good candidates for the warm component of the 
dark matter IIT4ll . Indeed, their only interactions with left-handed neutrinos are medi- 
ated by heavy bosons K and X (see Eq.©) via such couplings as -^{viiY"f^VLX^ 

and -^{v ji)^^^ V lK'^. When calculating the width for a kinematically allowed decay 
vn — !• v'j^vlv'j^, the expression F = G^m^^m^^ /192Ti^m\ satisfactorily sup- 
plies the order of magnitude. The life times r = F"^ ^ lO^^s, provided the fact that 
both bosons are in TeV mass region. This result is by far greater than the estimated 
age of the Universe ~ 4.3 x lO^^s according to L1.5J, so the right-handed neutrinos 
as WIMP particles have a safe behavior from the stability viewpoint. A deeper inves- 
tigation of the cosmological and astrophysical implications exceeds the scope of this 
letter and will be performed in a future work. 

5 Summary 

In this paper we worked out the canonical seesaw mechanism in the framework of 
electro-weak SU {A) l^U {1)y extension of the SM, by simply constructing dimension- 
five effective operators as suitable direct products among scalar quadruplets existing 
in the model. Then we made use of the same generic Yukawa matrix for Dirac and 
Majorana couplings in order to obtain the plausible mass spectrum for the left-handed 
neutrinos. It came out in the range 10^'^ — lQ~^eV . For this purpose, an intermediate 
scale between SM and GUT was set by the cut-off A ~ 1.4 x IQ^'^GeV. The right- 
handed partners develop masses in dependence on the high breaking scale V of the 
model. If the latter one ranges in 1 — \QTeV domain, then the right-handed neutrinos 
can be seen as plausible warm Dark Matter candidates with mass in IQ^^keV . 

The advantage of our method over other approaches is that it exploits in a natural 
way all the ingredients supplied by the 3-4-1 model itself, without enlarging the scalar 
sector of the model, without invoking any new particle in the lepton sector or a fine- 
tuning procedure. Of course, a more detailed analysis can further be performed based 
on different hypotheses regarding the Yukawa couplings that can be subject to a new 
and appropriate flavor symmetry. 
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